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Abstract 



We apply the closed-time path integral formalism to study the medium effects 
of non-equilibrium gluon matter. We derive the medium modified resummed 
gluon propagator to the one loop level in non-equilibrium in the covariant 
gauge. The gluon propagator we derive can be used to remove the infrared 
divergences in the secondary parton collisions to study thermalization of mini- 
jet parton plasma at RHIC and LHC. 
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I. INTRODUCTION 



Much effort has been given to detect a new form of matter known as quark-gluon plasma 
(QGP) for decades. The relativistic heavy-ion collider experiments at RHIC (Au-Au col- 
lisions at ^/s = 200 GeV) and LHC (Pb-Pb collisions at y/s = 5.5 TeV) will provide the 
best opportunity to study such a state of matter in the laboratory. Perturbative QCD 
estimates that the energy density of produced jets and minijets might be larger than 50 
and 1000 GeV/fm^ at RHIC and LHC [|l],0| which is much larger than the required energy 
density to produce quark-gluon plasma. As QGP lives for a very short time (several fermi) 
in a small volume (~ 100 fm^) a direct detection of this phase is not possible. Hence vari- 
ous indirect signatures are proposed for its detection. The prominent among them are: 1) 
J/\& suppression [Q, 2) strangeness enhancement 0], dilepton and direct photon production 
PJ^. However, many uncertainties exist which make it difficult to claim the existence of the 
quark-gluon plasma. The main uncertainty lies in the lack of an accurate determination of 
the space-time evolution of the quarks and gluons produced just after the collision of two 
nuclei at RHIC and LHC. While the quark-gluon plasma during the equilibrium stage is 
described by Bjorken's hydrodynamic evolution equations, it is much more difficult to de- 
termine the space-time evolution of partons in the pre-equilibrium stage. An accurate study 
of the pre-equilibrium stage will determine the equilibration time and initial conditions for 
hydrodynamic evolutions in the equilibrium stage. This study is also very important for the 
accurate determination of various signatures of the quark-gluon plasma. 

To describe the space-time evolution of quark-gluon plasma at RHIC and LHC one 
needs to know how the partons are formed in these high energy nuclear collisions. The hard 
parton (jets and minijets) production can be calculated by using pQCD. The pre-equilibrium 
evolution of these hard partons can be studied by solving relativistic transport equations with 
secondary collisions among these partons taken into account [l|,[^-[TO[| . However, soft parton 



production can not be computed within pQCD formalism. There are coherent effects for 



the soft partons and they may be described by formation of a classical chromofield |[Tl|--[17 . 
For simplicity, we will only consider the evolution of minijet plasma which can be studied 
by solving relativistic transport equation: 

p^d,fix,p) = Cix,p), (1) 

with secondary collision among the partons taken into account. In the above equation 

d^P2 r d^ps r d^p4 
(2^0M ' (27r)3pO J (2^2 
[f{x,ps)f{x,p,){l + f{x,p)){l + f{x,p2)) - f{x,p)f{x,p2){l + f{x,ps)){l + f{x,p,))] (2) 



^i^^P) = I fo \3 / fo \3 / fO \3 \^(PP^ P3Pi)?5^{P + P2-P3- Pa) 



is the collision term for a partonic scattering process pp2 — > PzPi- In the above expression p, 
P2i and p3. Pi are the four momentum of the partons before and after the collision. Different 
partonic scattering processes which can be considered are like gg gg, qq — >• qq, qq qq 
and gq — gq etc.. As the gluons are the dominant part of the minijet production we only 
consider gg gg process in the following. The squared matrix element for this process is 
given by: 

ms,u,t)\' = ^[3-^-^-^], (3) 
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where s = {p + P2f = {ps+PifJ = {p-p^f = {p2-P4f,u = {p-pif = {p2-p^f 
are the Mandelstam variables. For massless gluon they are related by 



t = -^[l~cosei u = —[l + cosei (4) 

where 9 is the center of mass scattering angle which goes from ^ ^ for identical partons 
in the final state. When one puts |M(s,-u, £)p from Eq. in Eq. (|^) one encounters 
divergence in the collision term at small angle {9 —>■ 0) or small momentum transfer t —>■ 0. 
This divergence is inevitable as long as we use the free propagator in vaccuum to evaulate 
the Feynman diagrams. However, this infrared divergence can be removed when one uses 
medium modified propagators instead of the vacuum propagator to evaluate the collision 
term in medium |18-pO|,^,pT|. The medium modified propagtaors have been obtained in the 



thermal field theory for the case of an equilibrium plasma. However, these finite temperature 
calculations are valid only in equilibrium where there is a static temperature such as in the 
case of a heat bath and the system is isotropic. For realistic situations in the high energy 
heavy-ion collisions at RHIC and LHC the partons formed at the initial time are in non- 
equilibrium and finite temperature QCD is not applicable at this stage. One has to use 
closed-time formalism to compute various quantities in non-equilibrium [p^-|25|. In this 



paper, using closed-time path integral formalism, we derive the medium modified resummed 
gluon propagator in non-equilibrium to the one loop order in the covariant gauge which is 
necessary to obtain finite collision term to study equilibration of expanding minijet plasma 
at RHIC and LHC. 

The paper is organized as follows. In section II we describe the formulation of closed-time 
path integral of SU(3) guage theory. In section HI we derive the resummed gluon propagator 
to one loop level. We summarize and conclude our main results in section IV. 



II. CLOSED-TIME PATH INTEGRAL FORMALISM IN GAUGE THEORY 

We consider SU(3) pure gauge theory which is QCD without quarks. In high energy 
heavy-ion collisions at RHIC and LHC most of the parton formed are gluons. Hence we 
concentrate on gluons only. However, extending close-time path integral formalism to quarks 
is straightforward. As the two nuclei travel almost at a speed of light at RHIC and LHC 
the system is dynamically evolving and many quantities have to be formulated in Boost 
invariant way pGf . For this expanding system of partons we work in the covariant gauge in 



this paper. The QCD action without quark is given by: 

S = l d'x[-\F'^^'^ F;, - ^^{C^f - CfpI (5) 
where the gluon field tensor: 

pat^u ^ Qi^j^au _ guJ^a^, ^ g f^^^A^" A'"' , (6) 

the ghost Lagrangian density: 

Cfp = [d^C]D^^[A]C (7) 
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and the gauge fixing term given by 



Qa 



(8) 



In Eq. (1^) ^ is the gauge fixing parameter in covariant gauge. In the closed-time path 
integral formalism the gauge field A, the ghost field C and the corresponding sources j, x 
are defined in both the time branches. To make the formulas simpler we denote the fields 
and the corresponding sources by following common notations: 



Q= {A.C.C), 



J = {j,x,x)- 



(9) 



Denoting Q and J+,J the fields and the sources on the upper and lower branch of 
the time path, the in-in generating functional becomes: 



Z[J+,J_,p] = j DQ^ <Q+Mp\Q',U 



> e 



(10) 



where S[Q^] = S[Q'^] — S*[Q ] with r = +, — . More explicitly, the above equation can be 
written as: 

Z[J+, J-,p] = J [dQ^]d[Q-] < g+,to| P \Q-,to > e'[^°[«l + ^-[Ql + ™ (11) 
where 

SoiQ] = 



E 



d^p A 



4 ,-A;{p)[G;l{p)rAt{-p) + C^{p)[S-\p)rC%-p)], (12) 



and 

TrJ ■ Q 



d'^p 



(27r)4 '2 



Mi-pHip) + c+{-p)x^{p) + xH-p)c^ip) + (+--)]• 



(13) 



In the Eq. (|T2|) G^y{p) and S{p) are gluon and ghost free propagators and in Eq. 

s^ntU = s.ntm-suQ-]. 

In covariant gauge the further complications arise because of the presence of the un- 
physical ghost fields. In equilibrium one can define the ghost distribution function (BE) 
and hence can deal with ghost fields in the medium even in covariant gauge. However, in 
non-equilibrium there is no easy procedure to obtain a ghost distribution function in the 
QCD medium at RHIC and LHC Quark and gluon distribution functions in non-equilibrium 
situations at RHIC and LHC can be obtained from minijets by using pQCD or via other 
methods 0-|n],|l],|3,|8l 



For this reason we work in the Landshoff and Rebhan scheme of 



frozen ghost formalism p9| , p0| , where the gauge theory in the covariant gauge is obtained by 
restricting the space of initial state \Q^,tQ > to the physical one. This means gluons with 
spatially transverse polarization will contribute to the trace. Hence the non-local Kernel 



K appearing in the path integral of the generating function |25| will couple only to the 
transverse component of the gauge field A. 

Now we consider a cylindrically symmetric expanding system of partons in 1©1 dimen- 
sion. For this purpose we introduce the fiow velocity of the medium 
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(cosh?7, 0, 0, sinh?7), 



(14) 



where 77 = | In ^ is the space-time rapidity and u^u^ = 1. We define the four symmetric 
tensors 0,11,11 



CM 



{u ■ pY — p"^ 

{u-p)pA I 



{u ■ pY — p^ 

1 

^2[{u-pY-P^] [ 

PfiPu 



U,. 



U, 



{u ■ p)Pf, 

p2 



_ {u ■ p)p^ 

p2 

Pu+ [Ui, 



{U-p)Py 

p2 



p^ 



and 



(15) 



which are required for a cyhndrically symmetric system in 1©1 expanding plasma at RHIC 
and LHC in the very early stage. 

Here T'^'" is transverse with respect to the flow-velocity but L'^'^ and D^'^ are mixtures 
of space-like and time-like components. These tensors satisfy the following transversality 
properties with respect to p^: 



P.T'^'ip) = P,L>^''{P) = 0, P^PuO^'ip) = 0. 

In addition to this, the above tensors satisfy the following properties: 

T ■ L = T -C = T ■ D = 0, T + L + D = l, 
T-T = T, L-L = L, C -C = \{L + D), 
TrC ■ L = TrC -0 = 0. 

Any symmetric tensor S'^" can be written in terms of the above four tensors: 

g^^u ^ a T^"" + b L^"" + c C^"" + d Df"" 

with 

^ -Tr C ■ S, d = Tr D ■ S. 



(16) 



(17) 



(18) 



a = -Tr T ■ S, b = Tr L ■ S, c 
2 



(19) 



In terms of this tensor basis the gluon propagator in the covariant gauge is given by: 



G^u{p)i 



iT^mpW + [G{P)] 



{9,. + (e - i)d,Ap)) [G{p)]7 - ^t,Ap) [G{p)]T;' ■ 



med 
ij 

VI 

ij 



i{L,,{p)+iD,,{p))[G{p)]\ 



where i, j = +, — . The forms of G{p)ij are 



[G{P)] 



vacuum 
ij 







2ni6{p^ 







'-Po) 



eipo) 



(20) 



(21) 



and 
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= -2mSip')f{p) i^Wy m = f{p) . eipo) + fi-p) ■ ei-po). (22) 

Here f{p) is the distribution function. Note that in this Landshoff-Rehban scheme the 
transverse component of the gauge propagator proportional to T'^'^ contains the medium 
effect which is gauge-parameter independent. As the initial density matrix contains only 
sum over transverse polarizations of gluons the ghost propagator is simply given by: 

S{ph = -I [G{p)]:;' , (23) 
because ghost fields do not couple to the kernel K in the generating function. 



III. RESUMMED GLUON PROPAGATOR AT ONE LOOP LEVEL IN 

NON-EQUILIBRIUM 



In the last section we derived the gauge field propagator in the medium for the free part 
of the action 5*0. In this section we will consider the interaction term of the action Sint- For 
such a situation one will have to consider the full propagator instead of the free propagator 
in the medium. This implies one has to solve the full Schwinger-Dyson equation which is 
practically impossible. Usually one has to truncate the series at one or two loop level of the 
self-energy. Such a truncation demands that the QCD coupling constant should not be too 
large. We consider here the evolution of the gluon-minijet plasma at RHIC and LHC where 
the average transverse momentum of the partons are found to be large at the earlry stage 
Therefore the coupling constant corresponding to such average transverse momentum is 
found to be small. In this situation we truncate the Schwinger-Dyson equation at one loop 
level of the self energy and consider the corresponding resummed two-point Green's function 
G. 

The resummed two-point Green's function G of gluon can be decomposed as: 



G,,{p) = -iT,,{p)G' (p) - iL^,{p)G^{p) - t^D,,ip)G^{p) 



(24) 



where G'^, G^, G^ correspond to T, L and D components respectively. The last part G^{p) 
is identical to the vacuum part ||3l| and hence we do not consider it any more. It can be 
mentioned that there are separate Dyson- Schwinger equation for different components which 
does not couple with each other. 

Dyson- Schwinger equations for different components can be written as 



G^'^ip) 



G^'Hp)].+T.[G^''iP) 

l,k 



n^'^(p) „ ■ a^'^^ip) 



Ik 



kj 



(25) 



Please remember that i,j, k,l = +,— and suppression of Lorentz and color indices in the 
above equation is understood. 

Instead of using the matrix form of this equation, we prefer to use the retarded, advanced 
and symmetric Green's functions. The retarded and advanced resummed Green's function 
are found to be: 



GVaIp) = G'r!a{p) + Gr^{p) ■ n^'^(p) ■ G^'^(p) 



yT,L 



tT,L . 



(26) 
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The straightforward solution of the above equation gives: 

yT,L , 



tT,L 



T.L , 



(27) 



where the self-energy contains the medium effects. Similar but more complicated equation 
is obtained for the resummed symmetric Green's function: 

+GYip) ■ ^YiP) ■ gY{p) + gY{p) ■ np^(p) ■ gY{p). 



(28) 



After some algebra it can be shown that 

Gl^\p) = [1 + 2f{p)] sgn{p,)[GY{p) - G^ip)] 
+ (np^(p) - (1 + 2f{p))sgn{p,)[lf/{p) - ^Ym) X G^ip) x GY{p). 



(29) 



Let us evaluate the various components of the self energy. For short hand notation we 
denote Gl for gluon loop, Ta for tapole diagram (see Fig. 1) and FP represents ghost loop 
diagram. Explicitly the expressions for various self energies are given by: 

,2 



n 



ab, 

Gl,fj,u 



(P) 



kl 



[G'-'mkl G^\p-q) 



n 



ah 

Ta;fj,i/ 



n 



ab 

FP;fj.u 



ip) 
(P) 



[9t^a{-p - q)i3 + gaf3{2q - p)^ + gf3f,{2p - q)c 
[9u^{p + q)s + g-tsip - 2g)z. + gus{q - 2p)^] . 
d^q 



kl 



kl 



kl 



-g NcSabSkl 

(k X l)tg^SabN, 



(2vr) 



■{gtixQucj - 9ax9^w)G'' {q)ki 



[G^'''mki[G^'''{p-<l)]kr<lM-p)u- (30) 



It has to be remembered that k,l{= (+,—)) are not contracted in the right hand side of 
the above equations. So it is obvious that ^T'a]^iu{p) f^i — 0' when k ^ I. In the following 
calculation we neglect the vacuum part and concentrate on the medium part of the self 
energy. Note that the distribution function is contained only in the medium part of the self 
energy. The divergence of the vacuum part of the self energy is absorbed in the redefinition 
of the bare quantities which is well known. Like the Green's function the self energy can be 
written as [|l|,||g: 

n^.(p) = T^.(p)n^(p) + V(p)n^(p). 



(31) 



Using Eq. ([TsD II's are decomposed as: 



n'^^^ (p) = KAp) 



n'^'^{p) = niM-L^''{p) 



-- R, A, S, 
K A, S. 



and 



(32) 



Total self energy of the medium part is the sum of gluon loop and tadpole contributions (see 
Fig. 1) because the ghost loop contribution is present in the vacuum sector in the frozen 
ghost formalism. The real and imaginary part of the total self energy can be written as: 



i?en^'^(p)„ = ReU'''^{p)Gl-a + n^'np)Ta;a 



R, A, S, 
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/mn'^'^(p), = /mn'^'^(p)Gi;a, a = R,A, S. 
One can check that the retarded and advanced self energy are related by: 

ReU^i%ip) = ^en^if^(p); ImU^i%{p) = -ImU^iAp) and 

U^iAp) = U^'Up)- (33) 

Furthermore we have ]l^'^.g{p)=0. Therefore we only need calculate RcHqiA Iml^Qi^, 
n^';.^ and IlTa;R which we compute in this paper for any non-equilibrium gluon distribution 
function fij)). Simplifying Eq. ( ^OD and using Eq. (|3^) the explicit expression of retarded 
self energy of the gluon loop diagram is found to be: 




FIG. 1. Gluon loop and Tadpole diagrams. 



where 
H^{q.p) 



+ 



{u ■ q){q ■ p){u ■ p) ^ {q-pf 



27r)3 {^p - qY isgn{pQ - qQ)e 

hp- q)^{{P - qf)H^'^{p - q,p) 
g2 + isgn{qo)e 

p^q-uf 



(34) 



{u ■ pY — p^ 



{u ■ pY — p^ {u ■ pY — p^ 



{p + qf 



{q-p){u-p) 



+ 



[q-pf 



+ 



P^ 



{u ■ q){{u ■ pY — p^) 2{u ■ qY{{u ■ pY — p^) 2{{u ■ pY — p'^) 



4^2 ^ g (g-p)(^-p) _ ^ {q-py 

{u -q) {u- qY 

- • Py + 2(g ■u)ip- u){q ■p) + {q- ufjp^ ~ 2{p ■ uY))) 

2(g ■ uYiiq ■ p) — p'^){p^ " {P ■ '^y) 
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(35) 



and 



H {q,p) = 7 ^ M • q) ^ 

[u ■ pY — p^ 

2{q ■ p){p ■ uY {q-pY{p-uY ■ m) 



{p + qy 



{q ■ u)p'^ {{p ■ uY — p^) {q ■ uYp'^Hp ■ uY — p'^) {{p ■ uY — p^) 



(m ■ g) (m ■ g)2 

. dp-uf{{q-pf-'^{q-u){p-u){q-p) + {q-uYp'^) . ^ 

^ {q ■ uYip"" - {p ■ uf) ■ ^ 

Using the 5 function we get: 



+ 



(27r)3 2|g1 (po -\q\Y -\q-p? + isgn{po - \q\)e 

'?o=-|'?l 1 



f{-q)-H^''^{q,p)\ 



(Po + 1^)^ - 1^- P\^ + isgn{po + \q\)e 

+ {q^p~q)- (37) 

Similar but more complicated calculation for symmetric self energy yields: 

- C^^^iv. / ((1,1 - p.f - 1,^- ,T) . ^^^'^^^^^'^^^'-'-^ ■ m • .1 

,2 ,^ i^(^-^)(g,p)|,o=-ki 



+ 5 [m + PoY - \q- P\') 2\T ■ ^^'"^ ■ ^^^^ 

where K^'^{q,p) are given by: 

r.Tf N p\p^ + 4:{q ■ u){{q ■ u) - {p ■ u))) ^ 3 

K {q,p) = — — -— _8g-u - lQ{p ■ u){q ■ u) 

4(g ■ u)^[q ■ u — p ■ uY\P — \P ■ uY) 

- I2{p^ ~ 2{p ■ uf) ■ {q ■ uf + A{'ip^{p ■ u) - A{p ■ ufYq ■ u) + S{p ■ uf + p'^ - Sp^{p ■ uf] (39) 
and 

K^q.p) = ^[8(? ■ - nP ■ n){q ■ uf 

4(g ■ uY\q ■ u — p ■ uYip — [p ■ uY) 

+ A(2{p-uf +p^){q-uf - V (P ■u){q-u)+ p'^] . (40) 

The expressions for the transverse and longitudinal part of the self energy coming from the 
Tadpole part (Eq. (|30|) and (|3^) ) are found to be: 

d^q f{q)u , {u-p){q-p) {q-pf 



(27r)3 2|gl ' 2(n • g)((u • pY - p'^) 2{u ■ qY{{u ■ pY - p'^) 
II , /(-^>ri , {u-p){q-p) {q-pf 



2((n-p)2-p2)Ji9o=|g| 2|g1 ^ ' 2{u-q){{u-pY -P^) 2{u ■ qY{{u ■ pY - p"^) 
2((n-p)2-p2)]l9o=-k1 (41) 



and 



^TaMP) = 9 SatN^J (^^[3 - 2 p2(^.^) + (^.p)2_p2(l " p2(^ . ^) ) ]I.O = k1 

+ 01 ^ ~ ^ — 27 \ ^ 7 ^2 2 2? '?o=-l'?l l^^i 

2\q\ p'^iu-q) [u-pY—p^ p^{u-q) ^ 

respectively. To summarize Eqs. (p7D, (p8D, (^Tj) and (^) contain all the expressions 
for different parts of the self energy in non-equilibrium which will be used in the medium 
modified propagator to obtain a finite collision term. In this paper we make hard momentum 
loop approximation in non-equilibrium which is simillar to hard thermal loop approximation 
in equilibrium In the hard momentum loop approximation the self energies are found 



to be gauge parameter independent (see Appendix). 

To obtain the collision term for the process gg — > gg in the medium we have to use 
the resummed Feynman propagator [G(p)++] which can be obtained from the resummed 
advanced, retarded and symmetric propagators via the relation: 

"^(P)]++ = \ [^s{p) + Ga{p) + GM] , (43) 

where 

GrAp) = -2 — — TTX^ — r^' (^^) 

- UrAp) ± ^s9n{Po)(^ 



and 



Gsip) = [1 + 2fip)] sgn{p,)[GR{p) - Ga{p)] 

+ {Iis{p) - (1 + 2f{f>))sgn{pomR{p) ' ^a^p)]) x Gr{p) x Ga{p). (45) 



Different parts of the self energy appearing in the above equations are given in Eqs. 
(^), (0) and (^2]). Finally, using the relation between various parts of the self energy (see 
Eq. (|33D ) we find from the above equations: 

^ . ^ P^ - Re^Ap) + \^s{p) . . 

{p^ - RenR{p)f + {ImnR{p)Y ^ > 

which is the required expression for the medium modified resummed Feynman gluon prop- 
agator in non-equilibrium at one loop level of the self energy. 

Our main purpose is to remove infrared divergence appearing in the small angle partonic 
scatterings which plays an important role in the production and equilibration of minijet 
plasma. For this purpose it is necessary to study the infrared behaviour of the self energies 
obtained in this paper for non-equilibrium situtations. We note that in the static limit 
(po = 0, IpI — > 0) one obtains n^(po = 0, |p| 0) = m|, (the Debye screening mass) and 
n^(po = 0, IpI — > 0) = TiYg (the magnetic screening mass). It is widely believed that while 
Debye screening mass is non-zero the magnetic screening mass is zero at one loop level of the 
self energy. For this reason one expects that while the electric field is screened, the magnetic 
field is not screened and one still has infrared divergence in the magnetic sector at one loop 
level. However, this is true for a system where the gluon distribution function is isotropic 
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(in momentum space) or in equilibrium. This is not true for non-isotropic gluon distribution 
function which is the case in the early stage of the heavy-ion collisions at RHIC and LHC 
The static limit result of the transverse part of the self energy derived in this paper is not 
zero for non-isotropic (in momentum space) gluon distribution function. The typical values 
of the magnetic screening masses at one loop level is found to be 257 MeV at RHIC and 
330 MeV at LHC by using non-equilibrium gluon-minijet distribution function at the 
initial time. Hence the medium modified gluon propagator derived in this paper is safe from 
infrared divergences both in electric and magnetic sector in non-equilibrium situations at 
RHIC and LHC. Therefore these propagators can be used to obtain finite collision terms to 
study equilibration of minijet plasma at RHIC and LHC. Note that the above values are 
obtained by using the minijet distribution function at the initial time. The time evolution of 
these screening masses have to be determined by solving relativistic transport equations with 
Bjorken's boost invariance picture taken into account. As these involves extensive numerical 
work (see |T[|) we will report the self-consistent space-time evolution study elsewhere. 



IV. CONCLUSIONS 

In this paper we have applied the closed-time path integral formalism to study medium 
effects of the minijet plasma in non-equilibrium. In particular we have derived the medium 
modified resummed gluon propagator to the one loop order of self energy which is necessary 
to obtain a finite collision term to study equilibration of parton plasma at RHIC and LHC. 
These medium modified propagators have been studied in more detail in finite temperature 
QCD formulations. However, at RHIC and LHC the parton momentum distributions at 
early stage are anisotropic and finite temperature QCD formulations can not be applied to 
these non-equilibrium situations. This is because a parton inside the nucleus (which travels 
almost at the speed of light at RHIC and LHC) carries mostly longitudinal momentum 
before an ultra relativistic nuclear collisions. After jets and minijets are formed and suffer 
secondary collisions the isotropy between longitudinal and transverse momentum may be 
achieved. According to Bjorken's proposal p6| many quantities are expected to be expressed 



in terms of boost invariant parameters. For this reason and to have a covariant formulation 
we have worked in the covariant gauge which is a suitable gauge for expanding plasma. We 
give the result of the resummed gluon propagator upto one loop level of self energy in non- 
equilibrium. Furthermore it is shown that these propagators are infrared divergence free, 
both in electric and magnetic sector. 

The medium modified resummed gluon propagator we derived in this paper will be used 
to obtain finite collision term for the gg — > gg scattering process in non-equilibrium to 
solve the relativistic kinetic equation (Eq. (||)) to study equilibration of minijet plasma at 
RHIC and LHC. In this way one does not have to put ad-hoc values for the momentum 
transfer cut-off which crucially changes all the properties and hence determination of all the 
signatures of the quark-gluon plasma. In future we hope to use our results of the resummed 
gluon propagators in the collision term to study production and equilibration of the minijet 
plasma at RHIC and LHC. As the solution of the relativistic transport equation involves 
extensive numerical work (see |^) we will report it elsewhere. 
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Appendix A: The explicit forms of self energies in Hard Loop Momentum 
Approximation: 

We use hard momentum loop (HML) approximation in non-equilibrium which is equiva- 
lent to hard thermal loop (HTL) approximation in equilibrium. In the hard momentum loop 
approximation the loop momentum q is harder than external momentum p [Q. In HML 
approximation we find that the whole expressions are independent of gauge-fix parameters 
^. In the hard momentum loop approximation we find: 

1 r <fq I ^ fiq) 



{u ■ pY - p^ J (27r)3 2|gl ^ - p.f - \q- p\^ 
■ [8(m ■ q){q- p){u ■ p) — 4(g ■ pY — 'ip^{q ■ 

+ „^ , ^^ l^iu -q^q- P){u ■ p) - 4(g ■ pf - Ap\q ■ uf 



and 



u ■ q [u ■ qY 

+ {q^p-q), (47) 



Imnl,Jp) ^^a^AT,^--!--, I ^^sgnip, - {{p^ - - Iq-p^) {|| 
[8{u ■ q){q ■ p){u ■ p) — 4(g ■ pY — 4,p^{q ■ uY 
- 2(p + Ap ■ <,) + 2(i^^£^ - 

+ ■wn(po + m (m + p»f - Vi-fi') 4if 

[8('U ■ q){q ■ p){u • p) — 4(g ■ pY — Ap^{q ■ uY 

u ■ q [u ■ qY 

+ iq^p-q), (48) 

Similarly for longitudinal components we obtain: 

ReUaiM ~ irSabNcj——, — / 77^1^7- ^o [8p {{u ■ q) 



2 [u ■ pY — p^ J (2vr)'^ 2|g| (po — 1^)^ — I? — P 

[u ■ qjp"^ [u ■ qYp 

^ J {2nY m iPo + \q]Y-\q-p\'^ ' P' ^ 



2 



+ {q^p-q), (49) 
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and 



sgnipo - \q\)6 {{\q\ - p,f - \q - pf) [Sp%u ■ q) - ^^^LlMlAf) 



p^ 

(p ■ qY{u-pf ^ ^ {q-pf{u-p f 
u ■ q)p^ {u ■ qYp^ 



{p-q)'^{u-pY Al'PYhj^'P) 



2 



+ (50) 

For the symmetric part of the self energy we find: 

ng^sb) = ^^^^^^'^ (^^.Iy_p2 I ((1^1 -Pof - \<i-p\") ■ m^o -u-qfrn ■ Hp- q) 

+ S {m +Por - \q-p\') ■ m-uo -u-qm-q) ■ hp- q). (51) 
It is interesting to see that in the Hard Loop Momentum Approximation we have n[^^^(p) ~ 

ng?5(p)- 
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